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We study the time correlation function of a density field in two-dimensional driven diffusive sys- 
tems within the framework of fluctuating hydrodynamics. It is found that the time correlation ex- 
hibits power-law behavior in an intermediate time regime in the case that the fluctuation-dissipation 
relation is violated and that the power-law exponent depends on the extent of this violation. We 
obtain this result by employing a renormalization group method to treat a logarithmic divergence 
in time. 
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I. INTRODUCTION 

The anomalous time correlation of hydrodynamic 
modes has been studied for a long period. For an equilib- 
rium fluid, it is understood that this anomaly arises from 
nonlinear mode coupling effects yj. By contrast, there 
is no systematic understanding of the time correlation 
in nonequilibrium steady states (NESSs) far from local 
equilibrium. In particular, it is not known how violation 
of the fluctuation-dissipation relation (FDR) influences 
the time correlation. 

As the simplest example realizing a NESS far from 
local equilibrium, we consider a two-dimensional driven 
diffusive system, in which a fluctuating density field is 
driven locally by an external force. Such a system can 
be realized in laboratory experiments Q. Perhaps the 
simplest model for a theoretical study of the long time 
behavior in the driven diffusive system is a stochastic dif- 
ferential equation consisting of terms representing a drift 
due to the external force, diffusion and random noise. 

The time correlation function for such a stochastic 
model has been calculated by employing mode coupling 
theory 3]. However, the model analyzed in Ref. |3J does 
not exhibit the long-range spatial correlation that is a 
generic feature of NESSs in driven diffusive systems of 
d > 2 dimensions. The reason that long-range correla- 
tion does not appear in that model is that violation of 
the FDR is not taken into account. Indeed, it is known 
that, in general, long-range correlation cannot exist when 
the FDR holds. By contrast, it has been found that the 
long-range correlation of driven diffusive systems can be 
described by a linear model with the violation of the FDR 
i. 

With the above considerations, in the present paper, 
we study a nonlinear model in which the FDR can be vi- 
olated. We demonstrate that as a result of this violation, 
the time correlation is qualitatively altered. Specifically, 
by employing a perturbative renormalization group (RG) 
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method that treats a logarithmic divergence in time |5|, 
we obtain an expression for the time correlation function. 
From this expression, we find that power-law behavior 
appears in the time correlation if and only if the FDR 
is violated and that the power-law exponent depends on 
the extent of the violation. 



II. MODEL 

We consider the time evolution of a fluctuating density 
field p(x,t) in a two-dimensional space, under the influ- 
ence of an external driving force in one direction, say the 
x\ direction, where x = (x\,X2). Note that we study 
NESSs in the high temperature regime, far from the crit- 
ical point. We now describe the model we study. First, 
the conserved quantity p obeys the continuity equation 

dp(x,t) A dJj(x,t) = 

dt ^ dx t 1 ' 

i—l 

We assume that the i-th component of the density cur- 
rent, Ji{x,t), is given by 

Ji(x,t) = -D i d i p(x,t) + S n J(p(x,t))+^(x,t). (2) 

Here, the functional form of J is such that, with p the 
average density, J(p) is the average current along the 
x\ direction in the steady state. We then approximate 
J(p(x,t)) in the form 

J(p(x, t)) » J(p) + c(p)Sp(x, t) + X(p)(5p(x, t))\ (3) 

where p(x,t) = p + Sp(x,t). The term ^(s, t) in (0 rep- 
resents a random current constituting zero mean Gaus- 
sian white noise, with 

(&(aj,t)&V,0) =2B t 5 tJ 5(x-x')8(t-t l ). (4) 

Note that because anisotropy in both the diffusion and 
noise intensity is expected to arise through effects of the 
external driving, the diffusion constant, Di, and the noise 
intensity, Bi, are assumed to be anisotropic, in general. 

Let us simplify the model given above. First, note 
that the first term on the right-hand side of (0) does 
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not contribute to the time evolution of the density, and 
the second term can be eliminated when we study density 
fluctuations in a frame moving with the velocity c given in 
©. To make this explicit, we define the density <j)(x, t) = 
5p{x + ce\t,t), where e\ is the unit vector in the x\ 
direction. Furthermore, introducing the parameters x 
and A, we rewrite B\ and Bi as 



Bi 



Ax(l-(-l) l A). 



(5) 



Thus, A corresponds to the extent of the violation of the 
FDR of the second kind @. Then, replacing Xi by \fD\xi, 
by y/xiPiDtfly and & by Vx^^M" 1 ^, we 
obtain the following dimensionless form of the equation 
for 6: 



d<p(x,t) 



E 



d 2 0(x,t) 0&(x,t) 



dxj 



dxi 



A 



(6) 



Here. 

(&(!c, *)&•(*', 0) = MiiQ- - (-l) l A)<5(a; - x')S(t - l'). 
and A is a dimensionless constant given by 

A = A(Di-D 2 )~ 1/4 X 1/2 - (8) 



(7) 



The renormalization group flow of (A, A) for the model 
© with is studied in Ref. 0- Also, the time correla- 
tion function has been calculated in the special cases that 
A = (using the mode coupling equation) [3| and A = 
P|. However, as far as we know, the time correlation 
function for the nonlinear model I© with the anisotropic 
noise intensity Q has never been investigated. 

In the analysis below, employing a perturbative ex- 
pansion with respect to A and A, we calculate the time 
correlation function C(k, t) defined by 



(2ir) 2 C{k,t)5(k + k') = U(k, 0)$(k',t) 



(9) 



Here and below, for an arbitrary function f(x, t), we de- 
fine f(k,t) by 



d 2 xe- lka 



f(x,t). 



(10) 



From the definition © and the symmetry of the steady 
state with respect to translation in time, the equality 



C(k,t) — C(k,—t) holds. Therefore, we consider only 
C(k,t) with t > 0. 



III. ANALYSIS 



First, we fix A and consider the expansion 
in A: 



of j>(k,t) 



4>(k,t) = ^ 0) (k,t) + X4>^(k,t) + X 2 ^ 2) (k,t) + --- . 

(11) 

Substituting Ijllfl into © with I|1U|) and extracting all 
terms proportional to A", we obtain a linear differen- 
tial equation for cjP 1 ^ containing all lower order and 
£i(k,t). Solving these differential equations under ini- 
tial conditions set at t = —00, we can iteratively derive 
expressions for (j>(°\ (j>( l \ ■ ■■ ■ We then substitute these 
results into ©• In this way, the correlation function 
C(k,t) is calculated in the form 

C(k,t) = C {0) (k,t) + XC {1) {k 1 t) + X 2 C {2) (k,t) + ■ ■ ■ . 

(12) 

It turns out that it is simplest to obtain the terms 
C (n ^(k,t) in the above expansion of C(k,t) by first de- 
riving the terms C^ n \k,uj) in the analogous expansion 
of C{k, u>), the Fourier transform with respect to time of 
C(k, t), and then taking the inverse Fourier transform of 
these. 

The lowest-order contribution to C(k,t) can be easily 
calculated as 



C(°>(fc,i) 



1 + A 



k 2 



k 2 



\k\ 2 



Ifert 



(13) 



Note that the spatial correlation function, obtained 
through the Fourier transformation of (k, 0) , exhibits 
power-law decay of the type 1/r 2 , unless A = 0. This 
illustrates the long-range correlation of driven diffusive 
systems. To this order, we find that there is no interest- 
ing behavior of the time dependence of C^(k, t), which 
merely exhibits an exponentially decaying form. 

The next contribution to C(fe, t) appears at second or- 
der in A. Through a straightforward calculation, we ob- 
tain 



6W(k,t) 



dt' 



d 2 k' Eti(i - e -ifen*-*'HfcTi«'i-ifc-fcTit'i 



(27T)* 



,'12 



• £j =1 (i-(-i^ 

\k\ 2 kl[kl 



^^-^w\ +]t - t,]+ w)-2W\ 



1 k 2 E-=i(l "(-IK AXfc; 



■j ^2 



Ifc- k'\ 2 



(14) 
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(In Sec. El we present a calculation method to obtain 
this result.) For this expression, we first perform the 
integration over \k | and then consider the t' integration. 
Next, we carry out the integration over the angle of k! . 
However, this procedure is complicated by the fact that 
divergences appear in the t' integration. As one example, 
(fc, t) includes the term 



- — kfte-W 2t tdt'-e^ 
^ 1 J V 



t'/2 



(15) 



where the contribution to the integral around t' = 
yields a logarithmic divergence. Physically, this diver- 
gence arises from the interaction between different modes 
during a very short time interval. However, the model 
we study is assumed to be appropriate only for describ- 
ing phenomena over time scales longer than a certain 
scale r m in driven diffusive systems |8(, and this diver- 
gence should not exist in the case that we study a model 
that correctly describes the phenomena with time scales 
shorter than r m . However, here, instead of studying a 
model in which the microscopic details of behavior on 
such short time scales are taken into account, we simply 
introduce a cut-off r m ; that is, the integration range of t' 
in ijbf|) is replaced by [— oo, — r ro ] U [r m , oo]. 

With the cut-off introduced, the term l|15|l can be reg- 
ularized as 



1 

47r' 



-kite 



\k\ 2 t 



log- 



t 



(16) 



Here, the first term exhibits a logarithmic divergence as 
t/i~ m — » oo, with fixed \k\ 2 t. Following similar proce- 
dures, we can separate all singular terms from C^ 2 \k, £), 
and in each case we obtain a term ~ \ogt/r m . 

Next, we expand (7 (2) (fc,i) in A. Then, to the first 
order, we obtain 

C(k,t) = C* (0) (fc,i) 



1 - (co(fe)A + Cl (A)k 2 t)\ 2 log — 

T m 

-A 2 C( 2 )(fc,<) + o(A 2 ,A), 



(17) 



where C' 2 ' (fe, t) represents the non-singular contribution 
to C^ 2 '(k,t), and co(fe) and Ci(A) are given by 



c (fe) 
C1 (A) 



k 2 k\ 



3k 2 



87r|fc| : 



A). 



(18) 
(19) 



(In Sec. [EJ we present a more detailed explanation of the 
derivation.) Note that the equal-time correlation C(k, 0) 
must be obtained as lim Tm ^o C(k, r m ), because the ex- 
pression 117fl is physically sound only for t > r m . 

The bare perturbation result (|17fl is reliable only for 
values of t for which logi/r m is of order unity. Now, 
employing the RG method demonstrated in Ref. we 



derive a form of C(k, t) reliable even for t ^> r m . First, we 
introduce a time scale tm which can be chosen arbitrarily 
and define a dimensionless parameter /i = tm/t u1 . Then, 
using 



log — = log — 

T m T M 



log 



(20) 



we rewrite (|17|l as 

C(k,t) = Z(n)C (0) (k,t) 



l-(c (k)A + Cl (A)k 2 t)\ 2 \og — 

tm 

+ \ 2 & 2 \k,t)+o(\ 2 ,A), 



(21) 



where we have introduced the renormalization constant 
Z(n). Here, the bare perturbation result l|17|) is equiva- 
lent to 121) with 

Z(fi) = l-(c (fe)A + ci(A)fc 2 i)A 2 log A1 + o(A 2 ,A). (22) 

Now, we regard (1221) as the bare perturbation result for 
Z(p) and calculate the improved perturbation result by 
using the fact that C(k, t) does not depend on tm- That 
is, differentiating (|21() with respect to tm, we obtain the 
equation 



dlogZ(n) 
dlogfi 



(co(fe)A + Cl (A)fc 2 i) A 2 + o(A 2 , A) = 0, 



which is referred to as the "renormalization group equa- 
tion" . Solving l|23|) under the condition 



we derive 



Z (fj,) = /i 



Z(/x = 0) = 1, (24) 

-(c (fe)A+ci(A)fe?t)A 2 +o(A 2 : A) ^5) 



which provides the improved result of l|22|l . Finally, sub- 
stituting (|25|l into l|21l) and setting tm = t (recall that 
tm is arbitrary), we obtain the expression 

/ , x -(c n (fc)A+ Cl (A)fe 2 t)A 2 +o(A 2 ,A) 

C(k,t) = C {a) (k,t) - 

+ \ 2 C^(k,t)+o{\ 2 ,A), (26) 
which may be reliable for all t > r m . 



IV. RESULTS AND REMARKS 

From the expression J2HJ, we have the following phys- 
ically interesting results. First, we note that there is a 
crossover time r c (fc) given by 



A|c (fe)|A 2 = |fe| 2 r c (fc). 



(27) 



We focus on the small wavenumber regime satisfying 
T c (fc) > r m . Then, for t satisfying t/\og(t/r m ) < r c (fe), 
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the correlation of density fluctuations takes the power- 
law form 



C(k,t) 



1 + A 



P - h 2 
Ifcl 2 



-c (k)A\ 2 



(28) 



It is important to note that this power-law regime 
appears only when A =/= 0, that is, only when the 
fluctuation-dissipation relation is violated [see JHJ)]. We 
believe that such power-law behavior can be observed in 
experiments. 

In addition to the above result, from l|26|l . we find 
that the decay rate of the correlation for t satisfying 
t/log(i/r m ) ^> T c (k) is expressed by 



t 



log C(M) - 1*1 + c 1 {A)k(X 2 \ogt. (29) 



This shows that the decay rate of the correlation increases 
slowly as a function of time. Such enhancement of the 
decay rate exists even in the case A = 0. A similar 
result was obtained from analysis of the mode coupling 
equation 0. 

The appearance of the singular term logi/r m in the 
bare perturbation result is the key to obtaining the 
power-law behavior of the correlation. A similar singular 
term was treated in Ref. |5( within the framework of the 
RG method to derive a solution representing anomalous 
diffusion for a deterministic nonlinear diffusion equation. 
There are several related works 0, in which such a 
divergence is treated in a similar way. 

The RG analysis given here should not be confused 
with the method to study the RG flow of system pa- 
rameters that occurs with the change of the wavenumbcr 
scale. For the model under consideration, this type of 
the RG flow of (A, A) is investigated in Ref. 7]. With 
that method, for example, the relevancy of the parame- 
ters can be studied, but explicit calculation of the time 
correlation is not possible. 

In conclusion, we calculated the time correlation func- 
tion 12tifl for the driven diffusive model © with Q. The 
expression we obtained indicates that a power-law regime 
appears in the time correlation function if the FDR is vio- 
lated. In addition to predicting this new type of physical 
phenomenon, our analysis provides an instructive exam- 
ple for the application of the perturbative RG method. 
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perturbative calculation of the time correlation function 
for the driven diffusive system. This work was sup- 
ported by a grant from the Ministry of Education, Sci- 
ence, Sports and Culture of Japan (No. 16540337). 



APPENDIX A: DERIVATION OF itTHl 

For an arbitrary function f(x,t), we define f(k,u>) as 

f(k,u) = J d 2 xdte- lut - lk x f{x,t). (Al) 
Then, the quantity C(k,u>) satisfies 

(2^) 3 <5(z + z')C(z) = (4>(z)4>(z>)) , (A2) 

where z = (k,uj). Here, the Fourier transformation of 
© yields 



(z) = G{z) 



with 



-y^jkj£i(z) - \iki{(j)o 4>){z) 



G{z) 



+ £i=i k2 ' 



, (A3) 



(A4) 



where (fog)(z) denotes the convolution of /(z) and g(z). 
From g3j, for $ n >(z),(n = 0,1,2,...), defined by |rfj) 
and IjAlfl . we obtain 



#°>(z)=G(z) 

. »=i 

(1) (z) =G(z) [-^i(0 (o) o ^ 0) )(z) 
4> (2) {z) =G(z) \-2ikx{ft 0) o<p))(z) 



(A5) 

(A6) 
(A7) 



We expand C(z) in the form 

C(z) = C (0) (z) + \C {1 \z) + \ 2 & 2) {z) + ... . (A8) 
The lowest order contribution of HASH is calculated as 



C<°>(z) = 2\G(z)\ 2 5> 2 (1 + (-lf-^A). (A9) 



i=i 



Using the inverse Fourier transformation in w, we obtain 
(Q3J. It can be easily checked C^(z) = 0, and C {2 \z) 
is expressed in the form 



C( 2 \z) = C?\z) + C^(z) + c£\z) 



(A10) 



where 
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c[ 2 \z) =8|G(z)| 2 fc 2 f d 3 z'\G(z z>)f J> ~ fc 2 (l - (-l) l A)|G(*')| 2 £ kf (1 - (-l^A), (All) 
Cg\z) =32|G(z)| 4 a,A ;i J2 ^(1 " ("I) 4 A) / dV|G(z - ^')| 2 E(% - k' 3 f{\ (-iy A)|G(z')l Vfei, (A12) 

i=l j=l 

G^) = - 32|G(^)| 4 |fc| 2 fc 1 X; *f (1 - (-1)*A) / d?z'\G{z - z'^j^ik, - k' 3 ) 2 (l (-iyA)\G(z')\*\k'\% 



(A13) 



Note that all the functions Ca {z), (a = I, H, HI), take 
the form 



C™{z) =F a (z)(h a ol a )(z), 



where F a (z),h a (z) and £ a (z) are determined from 
(IA11|) - (|A13I) . Using this form, we can express 
as 

£< 2) (M) = / dt'F a (k,t-t') 
d 2 k' 



APPENDIX B: DERIVATION OF (T7h 

(A14) ^ 6 ex P an< ^ C^{k,t) in the form 

G (2) (fc,i) = C^°\k,t) + AC^ 1} (k,t) + o(A). (Bl) 
Through a straightforward calculation, we obtain 



C^(k,t) = ^kf 



f dt g*|V/a_ i f* dt >l\k\H'/2 
Jo Jo *' 



(B2) 

^aA a (fc-fc,t)^a(fc,t). (A15) 

In order to calculate C^ 2,1 )(fe, t), we extract terms pro- 
portional to A from 114|) . The obtained expression be- 



Substituting this into (|A10fl . we obtain ijTIjl. 



comes 



/ ^/ 2 fr' r/r /2 -i-' 2 P - P 
(t-t') / « t "<1 "<2 i ft l 



(2tt) s 



|fc 



.'12 



(fc x - ^)e-[l fe - fe 'l 2+ l fe 'l 2 ] t ' 



. 2 j£_ / jVp-lfcl 2 !*-*'! /" d2fc/ fc i 2 Z fc 2 2 -f|fc-fc'l 2 +|fc'l 2 llt'l 
7 J (2^) 2 |fc'| 2 

■2— f ^'p-^I 2 !*-*'! / d2fc/ l" fc i 2 - k 2 , fc i - k V \ ( , _u ) -[\k-k'\ 
2 \k\zJ Me J (2^) 2 [ \k'\* + |fe|2 J ^ fc ^ e 



! +|fe'| 2 ]|*'l 



(B3) 



We first evaluate the Gauss integrals in |fc'| and perform Combining (|B2|) and (|B4(I with l|13|l . we finally obtain 
the t' integrals with picking up singular terms. Then, (|17|l with (|18fl and (|19|) . 

(j( 2:1 \k,t) is obtained as 



J_ -|fc| 2 t_^l_ ^ fc 2 3fcg \ , J_ 
~8n |fc| 2 ^ |fe| 2 y r m 

+ (non-singular term). (B4) 
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